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PACS numbers: 68.35.Fx, 05.40.-a, 02.50.Cw 

In Ref. [1| , we investigated a two-dimensional interface 
that grows in one octant of the cubic lattice. Exploit- 
ing limiting cases and symmetry, we conjectured two ba- 
sic nonlinear equations of motion for the interface speed. 
Combining these equations allows us to fit the interface 
speed along the (1,1,1) diagonal perfectly, but this fit- 
ting requires an unnatural choice of fitting parameters. 
Aesthetics suggests that one of our elemental equations, 


Zt 


Z X 1 Zy 1 


1 - 


1 


= R, 


(1) 


accurately describes corner 3d interface growth. While 
the prediction from ([1}, w = 0.125, for the interface speed 
in the (1,1,1) direction accurately matches our measured 
value w = 0.1261(2), small discrepancies persist. The 
comment Q studies the same interface growth rules start- 
ing from a fiat interface perpendicular to (1,1,1). This 
work finds w = 0.12606(2), consistent with our numerics, 
but also in slight disagreement with the solution to the 
conjectured exact equation. 

We recently found other independent equations that 
satisfy the required symmetries. One example is 
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for arbitrary n. Setting n = 1 + log 3 (8u>) perfectly 
matches numerics; for w = 0.12606, n = 1.0077. How- 
ever, the conjecture (JIJ is aesthetically more compelling. 

It seems coincidental that the beautifully symmetric 
growth equation ([T]) should differ from simulations by 
less than a percent. It is also unsatisfying to reproduce 
the numerics with high accuracy by using an equation 
such as (10) in our original Letter or @, which contain 


unnatural fitting parameters. We speculate that system- 
atic effects in the simulations may generate small dis- 
crepancies with the prediction of ([lj. Also note that 
little is known analytically about long-lived transients in 
(2 + l)-dimensional Kardar-Parisi-Zhang growth A 
recent numerical study [J] reveals a similarly stubborn 
approach to asymptotics in measurements of scaling ex- 
ponents for KPZ growth models. On a similar note, it 
is conceivable that differences in height correlation func- 
tions between the flat hypercube-stacking model [5] of 
the comment and the curved corner interface that we 
examine may generate slight differences in the interface 
speed w. Such differences in interfacial statistics between 
flat and curved geometries have been found rigorously in 
analogous (1 + l)-dimensional growth models @. 

While we agree with 0] that the numerics slightly de- 
viate from the predictions of (jTJ), it seems rash to reject 
this simple equation of motion in favor of an unnaturally 
complex one that minutely improves the accuracy for the 
interface speed. The outstanding challenge, of course, is 
to derive the correct equation for the interface motion. 
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